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Abstract. - Based on the Robertson theory the nonlinear dynamics of general Ising sys- 
tems coupled microscopically to bath systems is investigated leading to two complimentary 
approaches. Within the master equation approach microscopically founded transition rates are 
presented which essentially differ from the usual phenological rates. The second approach leads 
to coupled equations of motion for the local magnetizations and the exchange energy. Simple 
examples are discussed and the general results are applied to the Sherrington-Kirkpatrick spin 
glass model. 



Introduction. - Ising spin models are probably the simplest interacting many particle 
systems to which many equilibrium and non-equilibrium phenomena found in nature can be 
mapped. In many cases the analysis of these models exhibit the basic features and conse- 
quently lead to a first understanding of the reality. Therefore numerous investigations em- 
ploying Ising models exist in literature and are the subject of research activities. In addition 
to the questions in physics various problems arising in other fields (biology, social science, 
information theory, etc.) are treated on the basis of Ising models. 

Ising spin systems do not exhibit an intrinsic dynamics. Thus the kinetics of these models is 
exclusively caused by couplings to bath systems of the environment. Forty years ago Glauber 

suggested a phenomenological master equation for a one-dimensional chain of Ising spins 
with nearest neighbor interactions. Generalized to other spin arrangements and to other types 
of interactions (for a review see 1^) the approach of Glauber has been nearly exclusively used 
to describe the dynamics of the Ising spin systems. Note that the above comment also applies 
to the numerous Monte Carlo simulations performed on various models in the last decades |3] . 

The transition rates of these phenomenological equations are required to satisfy detailed 
balance which guaranties relaxation to equilibrium for non driven systems. But this require- 
ment is not sufficient to fix the rates completely. Therefore ad hoc assumptions have to be 
employed which lead to the claim EI that kinetic Ising models should not be regarded as 
faithful representations of processes occurring in real systems. 

It is microscopic quantum-statistical investigations which can remove the arbitrarily of 
the phenomenological equations. The tools for such a purpose are known since nearly half a 
century ago 0] and many general approaches exist. The application of these general theories 
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in Ising models, however, are rare. A first approacli of Heims |5j uses a very special and thus 
artificial spin-bath interaction. Further approaches, both based on special realizations of the 
heat bath, are presented by Martin [HI and recently by Park et al. [?]• 

All these microscopic approaches find transition rates different from A recent investi- 
gation for the linear response of a general Ising model based on the Mori theory results in 
the same conclusion. It is the aim of the present letter to calculate the transition rates under 
general aspects and to work out the generic behavior for the general non-linear case. Thus 
neither specific Ising models nor specific realizations of the baths are presumed. The Ising 
spins may even interact with more than one bath to include problems like those treated in [Sj . 
Calculations are performed on the bases of the Robertson theory TO. 11, , although equivalent 
result are obtained by employing ^ or the Nakajima-Zwanzig approach |11| . 

With rare exceptions like the original problem of Glauber [Tl, analytical solutions of the 
master equations are not known and therefore approximations have to be employed. The 
typical procedure involves the elimination of fast degrees of freedom in the master equation. 
This leads to a reduced description and equations of motion result for the slow or relevant 
variables. The latter equations, in many cases, can directly be found by the Robertson theory. 
Thus both, the master equation and the equation of motion of such reduced descriptions can 
simultaneously be derived which justifies the above choice of the method. 

The microscopic system:. - A system of N quantum spins s.^ with s = ^ is considered in 
the presence of time dependent external fields Hi — Hi (t) . The spins interact via an arbitrary 
Ising spin-spin interaction Jij{= Jji) and are described by the spin Hamiltonian 

= -2 ^ H,st + H"^" , = -2 51 Jvsts'j (1) 

i i,j 

where Ja = is presumed. Both quantum spin i operators with Ti — 1 and Ising spins 
Si{= 2sf) are used simultaneously. 

The assembly of spins interact via spin bath interactions B^iSi with a system consisting of 
Nb adiabatically isolated baths 7i„ (^). Thus the total Hamiltonian of the system is described 

by 

a Oil 

(2) 

During the dynamic evolution of the spin system thermal equilibrium is presumed for each 
individual bath system Ti,^ leading to the statistical operator 

i?^ - exp(- )/Tr exp(- /3„Hf ) . (3) 

These requirements on the bath systems can be realized by sufficient large bath heat capac- 
ities and by faster relaxation of the bath system. Thus tb ^ ts where tb and Tg are the 
characteristic relaxation times of the bath and the spin system respectively. In general all the 
temperatures (3a are different. Due to couplings to additional reservoirs the /3q, may even be 
time dependent. Thus the f3a{t) - together with Hi{t) - are preset quantities. Both the /3a (i) 
and the Hi{t) are presumed to have negligible variations on the fast time scale rg. (^) 

(^)This interaction represents the most general one spin coupling to the bath. Note that the contribution 
-B?s| has no effect and that the coupling of |^ or similar two spin couplings will lead to spatial correlations 
via H^^ . Such a behavior may be possible but is certainly not the generic case. 

(^)Eq.J5J is an approximation which requires this presumption for the Hi{t). 
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With these requirements there is no need for an exphcit specification of the bath Hamil- 
tonian and the bath operators Bai. As shown below, it is just the absorptive part Xaii^) 
of the dynamic bath susceptibihty 

/•oo 



which enters the calculation for values lots ~ 1- Expectation values with respect to are 
represented by {■ ■ ■)b and denotes the Liouvillian defined by L^A = [TC^ , A]. It is assumed 
that (Bai) B — holds, which can always be achieved by a renormalization of the terms of 
the Hamiltonian. Note that lots ~ 1 implies ujtb 1 and the first term of the low frequency 
expansion of Xai can be used. Apart from interesting exceptions (compare 7^) this yields in 
the generic case 

Xati^) ~ const LU for w < r^"^ (5) 

where the constant of proportionality may depend on (3a. 

Calculation:. - Let A'' {Si, ... Sn) be a fixed but arbitrary set of variables which are 
functions of the Ising spins Si . This set is assumed to give a sufficient dynamical description of 
some physical question and will be called observation level according to II . Without further 
specification of the observation level at this stage the general approach of Robertson 'lO'.'llJ 
is applied employing the usual perturbation treatment up to the second order in Ti.^^ and the 
standard Markovan approximation 13 . As all variables A^ commute with this leads to 

7 />OC 

^AA')t^- rfT([C/(t,t-r)H^^,[H^^,^^-]]), , C/(t,t-T)=exp(-zT[i^(0+i^]) 
at Jo 

(6) 

where the Liouvillians and are associated to and to , respectively. The time 
dependent expectation values (. . .)f are performed with the generalized canonical statistical 
operator 

Rt = Rf R^ with Rf = exp(l^)/Trs exp(V) and V = Y1 ^t^'' C^) 

fc 

where the time dependent Lagrange parameters are implicitly determined by 

{A'')t = TrsA''Rf . (8) 

Note that the system of eqs.(|H)l and ((SJ represent a closed set of nonlinear differential eqs. for 
the expectation values {A'^)t and the Lagrange parameters Xf. Thus together with the initial 
values for the {A'')t (or with the initial values for Xf) the time dependence of all quantities 
are completely determined by this set of equations. 

Analog to eq.lO can be brought in a more convenient form. Using again Xai{^) = 
Xl,{-uj) and x'Li^) = ± T (e^^-- - 1) (Bj, ^(L^ ± ^) i?± )b this leads to 

(^'>t - Ea^ (ra.{ tanh(/3„i7. + /3aX,) - S^}Af^)^ (9) 



dt 



with 



= Ea^ {Ta^{taMPo.H, + /S^X,) - tauh j^j ) ^ (10) 



^ x'Li2H, + 2X,) 
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For all functions F of the Ising spins is defined as 

F^^ = ]p:i,S,F{S^,...,SN) (12) 

which enter in the calculation from [sf, F] — ^2sfF[iy The operators of the internal field at 
sites i are 

3 

In the step leading to ea. ((Tn|l the relation [T^I5] 

(S,0)t = (Otanh(y[,]))t (14) 

was used which holds for any operator O of the Ising spins not involving the spin Si . 

Eqs.Q and pU|) represent in a compact form the most general results of this work. These 
results are valid for any observation level which has to be specified for concrete physical 
problems. As already pointed out, all slow variables of the specific problem have been included. 
One can also see that via ea. dlUH a second requirement on the set of variables results. Being, 
in principle, arbitrary the initial value of i?^o assumed to have the special form of ea. (|lU|l . 
By an extension of the observation level, however, this second requirement can always be 
satisfied for any given initial state. 

Here two special observation levels will be discussed. First, the observation level is analyzed 
which leads to the master equations and contains all functions F{Si, . . . , Sj^) of the N Ising 
spins. Obviously all of the above requirements are satisfied. The second observation level is 
spanned by all Ising spins 5*1, ... , Sn and by the interaction Ti.'^'^ . This choice implies that 
all thermal equilibrium states of the spin system are possible initial states. Consequently this 
observation level leads to a dynamical description of a standard thermodynamic system and 
is therefore denoted as standard system in the following. 

Master equation. - Let the |<t) be the common eigenvectors of the Si where a = 
{(Ti, • • -(Tjv} denotes the configuration of the spin system with the possible values ai = ±1. 
Then the observation level is spanned by all the 2^ projectors |cr)((T|. Introducing the no- 
tation o-W {cTi,--- • -ajv} , the relation {\a){a\)^^ = ct,/2 (|cr) (cr| + |crW) (cr(') |) 
holds. Eq.© immediately leads to the master equation 

^^p'^ =-Y^(wrp''-wr'%-''') (15) 

i 

for the occupation probabilities = (cr|_Rf |cr) . The transition rates Wf are given by 



Wr^J2%^-''^ tanh(/3„i7r)) = ^ "^^f^'^C! ,^ (16) 
^ 2 ^ 4[exp(2/3„(T,i7f)- Ij 



with 



^-- 4^^) H:^H. + Y.J.,a,. (17) 



The rates are consistent with the results of [7||S] but disagree with the usual ad hoc 
assumption of phenomenological kinetic Ising models which uses constant values of coeffi- 
cient C^i independent of the configuration cr. Thus for a faithful representation of processes 
occurring in real systems the full cr dependence of the C^j has in general be used. 
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Standard system. - This observation level will be characterized by V" = —^ex'H'^^ + 
which leads to = A; + Xex^i. Introducing the internal local-field probability 
functions 

= {Six - X.)U = (18) 

and employing again the relation (|10() the expectation values of the magnetizations — {Si)t 
and the exchange energy [/"^^ — {'H'^^)t can be written as (compare |14ll5|') 



r..=ydx..(.) tanh(A.+A..a:) , C/^ ^ -i E J . -.(x) tanh(A. + A...) (19) 
and the cqs.© leads to 

=-E„/dxz..(x)^^^«±H^[,|^HM^^ (20) 

r tanh(A^+Aecca;) 



=E.Jdxxwdx V-''f^'^^ [ Zt^^^^^^ . (21) 



The set of nonlinear equations 1)1912111 for the dynamic variables rrii and W^^ together with 
the associated time dependent Lagrange parameters Xi and Xex is closed and is expected to 
give an adequate description for many physical systems and situations. 

Note that the internal local- field distribution functions Wi{x) govern the set of nonlinear 
equations of motion. Indeed all the terms of these set of equations can explicitly be calculated 
from the knowledge of the Wi{x). These findings generalize the previous results, that the Wi{x) 
determine the statics jl4| and the dynamical linear response j8j of an arbitrary Ising model. In 
that previous work the distributions Wi (x) have explicitly been calculated for various specific 
physical models including both ferromagnetic and disordered systems. Consequently for all 
those models the nonlinear equation of motion can easily be obtained from the eas. p9l21(l . 



Application to the SK spin glass. - The latter procedure will be illustrated for infinite 
range spin glass model of Sherrington and Kirkpatrick ^JEl- In the SK model the bonds 
Jij are independent random variables with zero means and standard deviations N^^. This 
scaling fixes the spin glass temperature to T = 1. Employing the modified TAP approach ^H] 

^SK(. „ 1 cosh(A. + Ae.x) ix-Hrr + {X,,A)\ 

' ^ ^" V2^cosh(A. + Ae.i?r) ^ 2A ^ 

holds IE] with 

\ /_■ _■ m ■ — A _ _. m - A Q nrl f n A — 

N ^l + r^XUl-m^) 



nr - E J^,m, - Ae.m.A and with A = 1 ^ n - ^2^2 ' (23) 



where F is determined from 



r = for 1-^ ^(1-^2)2 >o (24) 

i 



With the eas. (|22l25|l the equations of motion p9l21|l are explicit provided the system of baths 
is specified. For the case of one singular bath and identical spin bath couplings the sums over 
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a drop out in eas. H2UI21|) and only one bath temperature remains. The resulting equations 
are expected to describe the complete non linear dynamics of the SK spin glass including 
temperature quenches, memory and ageing effects |17| . 

A complete discussion of these interesting effects is beyond the scope of this work and is a 
subject of further research. In the following, I focus on the discussion of the special situations 
when the exchange system is in equilibrium with the bath, implying Xex = 13b- This situation 
can physically be realized by allowing only changes of the fields Hi. Setting Aex = (3b = (3 
ea. (|20|l yields Glauber type equations of motion 

m, = -7(/3, H^", A) [m, - tanh/3iJf^^] , H^" ^ + ^yJr.mj - m,/3A (26) 

with the relaxation rates 

7(/3, H, A) = / dz cxp <^ - - — } ^J— > (27) 

^ ' V2^ J I 2 i 2 sinh/3(i/ + \/Az) ^ ^ 

where A has to be determined from ea. H23l25|l (with X^^ replaced by /3). On a phenomeno- 
logical basis - with 7 replaced by a constant - equations of motions of this kind have been 
used since the early days in spin glass research Moreover the author himself had used 
such equations to find numerically the solutions of the static TAP equations jl8[l2U| . Note 
that ea. H27() implies relaxation to the static TAP equations independent of the values of 
7(/3, H^^^ , A). Thus the phenomenological equations can be used instead of the microscopic 
results to calculate static properties and these results justify the procedure used in (181120) to 
calculate the solutions of the TAP equations. 

For all dynamic quantities, however, the full form of ^{j3,H^^^ has to be used. First 
examples are the dynamic susceptibility and the relaxation function as presented already 
in 1201 ■ Note that a result of [20] , the frequency-dependent shift of the cusp temperature of 
the real part of the susceptibility, is basically a consequence of 7(/3, H^^^ , A). In this context 
the general result '11' should be recalled that the linear approximation of the Robertson theory 
leads the Mori theory. Consequently all the linear response theory results of can also be 
obtained from the of eas. ((26|l and H27|l by linearizing around the thermal equilibrium which 
can easily be checked. 

Discussion. - The consequences of the results of this letter are further illustrated on 
simple physical situations. First one singular spin coupled to a bath with temperature (3 is 
considered. This implies no spin spin interaction at all. Both the master equation and the 
standard system approach of this letter leads to the equation of motion m = — 7i(m— tanh/?_ff) 
for the magnetization m with 71 = ixB(2^^) coth/3i7 . For a static field the magnetization 
m relaxes to the equilibrium value tanh/3ff with the well known (see e.g. ^I]) longitudinal 
relaxation rate 71. Note the latter result is not found if instead of ea. H17|) the usual ad hoc 
assumptions of the kinetic Ising model are used. 

Next the original system of Glauber, an infinite chain of Ising spins with next neighbor 
ferromagnetic interactions is considered which is coupled to one bath with a temperature (3. 
The transition rates ((17)1 satisfy the detailed balance relation and therefore relaxation to the 
thermal equilibrium results within the master equation approach. The same consequences 
apply for the standard system approach. Indeed the equilibrium values of the Lagrange 
parameter Ai = Hi and Xex = (3 describe a fixed point of the eas. l)20() and l(21() . It is obvious 
that nether the transition rates 1(17(1 nor the details of distribution functions Wi enter in 
thermal expectation values. Note that these arguments apply to all undriven Ising system 
which are coupled to one singular bath. Thus all this systems relax to thermal equilibrium 
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taking in account that the fixed points are always stable. This can easily be checked by a 
linear stability analysis. 

For situations, however, where two or more bathes with different temperatures Pa are 
present the values of stationary solutions depend on the transition rates or on the details of 
the distribution functions Wi{x). Thus for example the non equilibrium stationary states of 
Schmittmann and Schmiiser 9 can not be accepted as a realistic solution for their interesting 
model of a one dimensional Ising chain coupled to two bath systems as the analysis is based 
on the phenomenological rates and not on the physical rates. 

Finally it is pointed out that the results of the present work will always effect the details of 
the dynamic evolution apart form situations where the hight temperature approximation can 
be used. In general more complicated expressions are found compared to the phenomenological 
results. This even applies to the original model of Glauber. For this kinetic Ising model the 
equations of motion for the time correlation function partly separate which enables the explicit 
calculation of the complete dynamics. Such a simplifying separation does not occur for the 
microscopical master equation. Thus the solution is not known and is hard to find as already 
pointed out in 1 . 
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